Lorenz number in high T c superconductors: evidence for bipolarons 
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Strong electron-phonon interaction in the cuprates has gathered support over the last decade in 
a number of experiments. While phonons remain almost unrenormalised, electrons are transformed 
into itinerent bipolarons and thermally excited polarons when the electron-phonon interaction is 
strong. We calculate the Lorenz number of the system to show that the Wiedemann-Franz law 
breaks down because of the interference of polaron and bipolaron contributions in the heat flow. 
The model fits numerically the experimental Hall Lorenz number, which provides a direct evidence 
for bipolarons in the cuprates. 



The discovery of high-temperature superconductors 
jl], ||] has broken constraints on the maximum T c pre- 
dicted by the conventional theory of low-temperature 
superconducting metals and alloys. Understanding the 
pairing mechanism of carriers and the nature of the nor- 
mal state in the cuprates and other novel superconduc- 
tors has been a challenging problem of the Condensed 
Matter Physics. A number of theoretical models have 
been proposed, which rely on different non-phononic 
mechanisms of pairing (see, for example ||, 0]). On the 
other hand, over the last decade, increasing evidence for 
the electron-phonon interaction has been provided by iso- 
tope effect measurements [H , infrared and ther- 
mal conductivity ||, neutron scattering |10| , and more 
recently by ARPES (ill p). 



To account for the high values of T c in the cuprates, 
one has to consider electron-phonon (e-ph) interactions, 
which are larger than those used in the intermediate cou- 
pling theory of superconductivity jl3) . Regardless of the 
adiabatic ratio, the Migdal-Eliashberg theory of super- 
conductivity and the Fermi-liquid have been shown to 
breakdown at the e-ph coupling constant Awl . The 
many-electron system collapses into the small (bi)polaron 
regime at A > 1 with well separated vibration and charge- 
carrier degrees of freedom. Although it might have been 
thought that these carriers would have a mass too large to 
be mobile, the inclusion of the on-site Coulomb repulsion 
and the poor screening of the long-range e-ph interaction 
do lead to mobile intersite bipolarons [^5, 16 1. Above T c 
the Bose gas of these bipolarons is non-degenerate and 
below T c their phase coherence sets in and superfluidity 
of the doubly-charged 2e bosons can occur. In this pic- 
ture, the thermally excited single polarons co-exist with 
the Bose gas. 

There is much evidence for the crossover regime at T* 
and normal state charge and spin gaps in the cuprates 
p7| . These energy gaps could be understood as being 
half of the binding energy A and the singlet-triplet gap 
of preformed bipolarons, respectively |Q. Many other 
experimental observations were explained using the bipo- 
laron model [|l9|. These include the Hall ratio, the Hall 
angle, ab and c-axis resistivities, magnetic susceptibility, 



and angle-resolved photoemission. The bipolaron model 
provides parameter-free fits of critical temperatures, up- 
per critical fields, explains a remarkable increase of the 
quasiparticle lifetime below T c |2l[| , and the symmetry of 
the order parameter p2| in many cuprates. Further ev- 
idence for bipolarons comes from a parameter-free es- 
timate of the renormalized Fermi-energy 1 23 , which 



yields a value well below lOOmeV. It is so small that 
pairing is certainly individual in most cuprates, i.e. the 
bipolaron size is smaller than the inter-carrier distance. 
This is the case in a (quasi) two-dimensional system, if 



e F < 2ttA. 



(1) 



The normal-state pseudogap, experimentally measured 
in many cuprates, was found as large as A/2 > 50meV 
p7[ , so that Eq.(l) is well satisfied in underdoped and 
probably also in optimally doped cuprates. One should 
notice that a coherence length in the charged Bose gas 
is not the size of a boson. It depends on the interparti- 
cle distance and the mean-free path, pl|, and might be 
as large as in the BCS superconductors. Hence, it is in- 
correct to apply the ratio of the coherence length to the 
inter-carrier distance as a criterion of the BCS-Bose liq- 
uid crossover. The criterion of real-space pairing is given 
by Eq.(l). 

Direct evidence for the existence of charged 2e Bose 
liquid in the normal state cuprate materials is highly 
desirable. In 1993 Mott and Alexandrov|f24[ discussed 
the thermal conductivity k; the contribution from the 
carriers given by the Wiedemann-Franz ratio depends 
strongly on the elementary charge as ~ (e*)~ 2 and should 
be significantly suppressed in the case of e* = 2e com- 
pared with the Fermi-liquid contribution. As a result, 
the Lorenz number, L = (e/ks) 2 K e /(Ta~) differs signifi- 
cantly from the Sommerfeld value L e = 7r 2 /3 of the stan- 
dard Fermi-liquid theory, if carriers are bipolarons. Here 
K e , a, and e are the electronic thermal conductivity, the 
electrical conductivity, and the elementary charge, re- 
spectively. Ref. [53 predicted a very low Lorenz number 
Lb for bipolarons, L b = 6L e /(47r 2 ) ks 0.15L e , due to the 
double charge of carriers, and also due to their nearly 
classical distribution function above T r . 
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Unfortunately, the extraction of the electron thermal 
conductivity has proven difficult since both the electron 
term, K e and the phonon term, n p h are comparable to 
each other in the cuprates. Some experiments have at- 
temped to get around this problem in a variety of meth- 
ods |^5|, |26|, |27j. In particular, Takenaka et al.|25|] found 
that K e is constant or weakly T-dependent in the nor- 
mal state of YBaiCu-sO^+x. This approximately T- 
independent n e therefore implies the violation of the 
Wiedemann-Franz law (since the resistivity is found to 
be a non-linear function of temperature) in the under- 
doped region. The breakdown of the Wiedemann-Franz 
law has been seen also in other cuprates J2^, |2!| . 

More recently a new way to determine the Lorenz num- 
ber has been realised by Zhang et al.[^o), based on the 
thermal Hall conductivity. The thermal Hall effect al- 
lowed for an efficient way to separate the phonon heat 
current even when it is dominant. As a result, the 
"Hall" Lorenz number, Lh = (e/kB) 2 K xy /(Ta xy ), has 
been directly measured in YBa2C'u30e.95 because trans- 
verse thermal n xy and electrical a xy conductivities in- 
volve only the electrons. Remarkably, the measured value 
of L xy just above T c is about the same as predicted by 
the bipolaron model, L xy w 0.15L e . However, the exper- 
imental L xy showed a strong temperature dependence, 
which violates the Wiedemann-Franz law. This experi- 
mental observation is hard to explain in the framework 
of any Fermi-liquid model. 

In this letter we propose a theory of the Lorenz num- 
ber in the cuprates explaining the experimental results 
by Zhang et al |3(|. Our particular interest lies on 
the conclusions that the Wiedemann-Franz law is vio- 
lated in the cuprates in the temperature range below the 
crossover temperature T* . Here we demonstrate that the 
Wiedemann-Franz law breaks down because of the in- 
terference of polaron and bipolaron contributions to the 
heat transport. When thermally excited polarons are 
included, the bipolaron model explains the violation of 
the Wiedemann-Franz law in the cuprates and the Hall 
Lorenz number as seen in the experiment. 

Thermally excited phonons and (bi)polarons are well 
decoupled in the strong-coupling regime of the electron- 
phonon interaction]!^], so that the conventional Boltz- 
mann equation for renormalised carries can be applied. 
We make use of the r— approximation |nj in an electric 
field E = V</>, a temperature gradient VT, and in a mag- 
netic field B _L E, VT. The bipolaron and single-polaron 
non-equilibrium distributions are found as 



/(k) = fo(E) 



dfo fv 



6n x F}, 



(2) 



k 2 /{2m p ) and 9 = Q p — eBT p /m p for thermally ex- 
cited polarons. Here m b ^ p are the bipolaron and polaron 
masses of too-dimensional carriers, y = exp(/i/T), /i is 
the chemical potential, h = c = ks — 1, and n = H/B is 
a unit vector in the direction of the magnetic field. Eq.(2) 
is used to calculate the electrical and thermal currents in- 
duced by the applied thermal and potential gradients as 



j a = a a/3 V/3(p - 2e4>) + bapVpT, 



(3) 



(4) 



Eq. (3) defines the current with the polaronic conductivity 
(j p = e 2 T p n p /m p , where the kinetic coefificents are given 
by 



1 

Ye 



a xx = a yy = — a p (l + L4), 



(5) 
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cr p (e p + L49 b ), 



J xx — w yy 



A - n 



e ' ' 2T 
2A(T b - fi/T)}, 

b yx — —b xy — —— [Q p ^T p H — 

2AQ b {T b - n/T)] 

Eq. (4) defines the heat flow with the coefficients given by 



2A(TT b + 2e<t>)}, 



c xx = c yy = ^[TT p + A/2 + eqb + (6) 



c yx = -c xy = -^[Q p (TTp + A/2 + eq 



2AQ b {TT b + 2e<p)], 



d xx = d yy = -|{T7 P + r p (A - fi/2 + eq 
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2T 



A[T lb + r fc (2e0 - n) - 2e<j>fi/T\}, 
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A - u, 
) - 

2T J 



Here 



where v =dE/dk, F = (E - /i)VT/T + V(/z - 2eqb) and 
fo(E) = [y~ 1 exp(E/T) — 1]~ for bipolarons with the an( j 
energy E = k 2 /(2m b ), and the Hall angle O = <d b = 
2eBn/m b , and F = (E + A/2-^/2)VT/T+V(^/2-e0) 
and f (E) = {<r 1/2 exp[(£ + A/2)/T] + l}" 1 , E = 



AQ b [T lb + r 6 (2e<£ -n)- 2e^/T)}} 

r _ J™dEE 2 df /dE _ 2$(z,2,l) 
Tj^dEEdf /dE $(2,1,1) 

_ / °° dEE 3 df /dE _ 6$(z,3,l) 
' ~ T 2 J 00 dEEdfa/dE ~ 1, 1) 



3 



are numerical coefficients, expressed in terms of the Lerch 
transcendent $(z, s,a) = J2'kLo zk /( a + k) s with z = y 
in T b , 7 b and z — — y 1 / 2 exp[— A/(2T)] in r p , 7 p , and 
A = m p T b n b I (m^TpUp) is the ratio of the bipolaron and 
polaron contributions to the transport, which strongly 
depends on the temperature. For simplictiy we neglect 
the spin gap, which is small in the optimally doped 
cuprates Il7| . Then the bipolaron singlet and triplet 
states are nearly degenerate, so that the bipolaron and 
polaron densities are expressed as 



rib = 1 hi(l - y)\, 



y 1/2 exp ( - A 



(7) 



(8) 



Using the kinetic coefficients Eqs. (5) and (6) we obtain 



P 
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a p {l + AA)' 
1 + AAQb/Sp 
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L p + AAL b A [2T p - T b + A/T} 2 



1 + 4A (1 + 4 A) 2 

L p + 4AL b e b /e p 

1 + 4Ae b /e p 

A(4A + & b /&p) [2T P -T b + A/T} 2 
(l + 4A) 2 (l + 4A9 b /e p ) 



(9) 
(10) 

(11) 
(12) 



for the in-plane resistivity, the Hall ratio, the Lorenz 
number and the Hall Lorenz number, respectively, where 
L p — (7 P — T 2 ) and L b — (7b — r 2 )/4 are the polaron and 
bipolaron Lorenz numbers. In the limit of a pure pola- 
ronic system (i.e., A = 0) the Lorenz numbers, Eqs. (11), 
(12) are L = Lh = L p . In the opposite limit of a pure 
bipolaronic system (i.e. A = oo) we obtain a reduced 
Lorenz number [p4| L = Lh = L b . However, in gen- 
eral our equations (11) and (12) yield the temperature 
dependent Lorenz numbers that differ significantly from 
both limits. The main difference originates in the second 
terms in the right hand side of Eqs. (11) and (12), which 
describe an interference of polaron and bipolaron contri- 
butions in the heat flow. In the low-temperature regime, 
T <C A, this contribution is exponentially small because 
the number of unpaired polarons is small. However, it is 
enhanced by the factor (A/T) 2 , and becomes important 
in the intermediate temperature range T c < T < T* . 
The contribution appears as a result of the recombina- 
tion of a pair of polarons into the bipolaronic bound state 
at the cold end of the sample, which is reminicent to the 
contribution of the electron-hole pairs to the heat flow in 
semiconductors |3l| . These terms are mainly responsible 
for the breakdown of the Wiedemann-Franz law in the 
bipolaronic system. 



It has been shown that the bipolaron model fits nicely 
the temperature dependencies of the in-plane |S2| and 
out-of-plane|3j| resistivities and the Hall ratio in the 
cuprates. Here we show that it also fits the Hall Lorenz 
number measured by Zhang et al.|30 . To reduce the 



number of fitting parameters we take the charge pseu- 
dogap A/2 = 600X, as found by Mihailovic et al.(f| 
for nearly optimally doped Y BaiCu-sO^+x in their sys- 
tematic analysis of charge and spin spectroscopies. Ac- 
cording to Ref. p(J the main scattering channel above 
T c is due to the particle-particle collisions with the re- 
laxation time Th.p oc 1/T 2 . The chemical potential is 
pinned near the mobility edge, so that y i=s 0.6 in a wide 
temperature range, if the number of localised states in 
the random potential is about the same as the number of 
bipolarons|2(J. This is the case in Y Ba^CuzO^x, where 
every excess oxygen ion x can localise the bipolaron. As a 
result, there is only one fitting parameter in Lh, Eq.(12) 
which is the ratio of the bipolaron and polaron Hall an- 
gles <d b /Qp. The model well fits the experiment, Fig.l, 
with a reasonable value of Q b /Q p = 0.44. It also quanti- 
tatively reproduces the (quasi)linear in-plane resistivity 
and the inverse Hall ratio, as observed in the cuprates 
(upper inset in Fig.l and Ref. §3, ||). 




T(K) 

FIG. 1: The experimental Hall Lorenz number j3^] in 
Y Ba2CuiO%.^ fitted by the bipolaron theory, Eqs. (11, 12), 
with Bb/0p = 0.44. The upper inset shows the linear in- 
plane resistivity and the Hall ratio. The lower inset shows 
the ratio of the Hall Lorenz number to the Lorenz number as 
a function of temperature. 

Lower inset in Fig.l shows a slightly lower Lorenz num- 
ber compared with the Hall Lorenz number produced by 
our calculations. Because the thermal Hall conductivity 
directly measures the Lorenz number in the framework 
of our model, it can be used to measure the lattice con- 
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tribution to the heat flow as well. When we subtract the 
electronic contribution determined by using the Lorenz 
number, the lattice contribution to the diagonal heat flow 
appears to be much higher than it is anticipated in the 
framework of any Fermi-liquid model. 

We notice that some recent measurements p4j on 
TI2BCL2CUO6+S suggest that the Wiedemann-Franz law 
holds perfectly well in the overdoped region and there- 
fore conclude that the Fermi- liquid prevails at this doping 
range. Alexandrov and Mott |l8j, suggested that there 
might be a crossover from the Bose-Einstein condensa- 
tion to a BCS-like polaronic superconductivity across 
the phase diagram. Thus Proust et al's results are still 
compatible with the (bi)polaron picture. If the Fermi 
liquid does exist at overdoping then it is likely that the 
heavy doping causes an " overcrowding effect" where the 
polarons find it difficult to form bipolarons due to the 
larger number of competing holes [fl8||. 

We conclude that by the necessary inclusion of ther- 
mally excited polarons as the temperature rises, the bipo- 
laron model predicts the Lorenz number very close to ex- 
periment in underdoped and optimally doped cuprates. 
Our consideration leads to good fits for the experimental 
Hall Lorenz number, Hall ratio, and the in-plane resistiv- 
ity. The interference of the polaron and bipolaron contri- 
butions to the energy flow breaks down the Wiedemann- 
Franz law and results in the unusual temperature depen- 
dence of the Lorenz number. 
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